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ABSTRACT 


An analytical investigation was conducted to determine the thermal degradation of 
high-performance multilayer insulation on a liquid hydrogen tank in the vicinity of a pene- 
tration. Results were obtained with and without the use of a thermal buffer zone between 
the shields and the penetration, for various values of thermal contact resistance. Input 
variables to the program included radiation source temperature, penetration diameter, 
number of shields, shield emissivity, shield thermal conductivity, shield spacing, buffer 
zone dimensions, and buffer zone thermal conductivity. The results of the program are 
presented in graphical form, as well as in equations representing curve fits of the ana- 
lytical data. 
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SUMMARY 

An analytical investigation was conducted to determine the direct lateral heat flow to 
a penetration, the net increase in heat flow into a tank, and the radius of the thermally 
degraded area created by the presence of a cylindrical pipe or tank support penetrating 
the multilayer insulation on a liquid hydrogen tank. The multilayer insulation system 
simulated by the analysis used nonconducting spacers between consecutive shields of 
1/4-mil (0.00064-cm) double aluminized polyester film. The analytical model considered 
the multilayer insulation butted directly to the penetration, as well as butted against a 
thermal buffer zone placed between the shields and the penetration. 

The resulting equations from the analysis were programmed for solution by a com- 
puter. Input variables to the program included radiation source temperature, penetra- 
tion diameter, number of shields, shield emissivity, shield thermal conductivity, shield 
spacing, buffer zone dimensions, and buffer zone thermal conductivity. Various values 
of thermal contact resistance between the shields and the penetration or buffer zone were 
generated to make the analytical results applicable to the real case. 

The results of the program are presented in graphical form to demonstrate the effect 
of any given variable on the thermal degradation of the insulation system. Specific values 
of the heat flows into the tank, for a particular penetration and insulation system, can be 
determined for any value of thermal contact resistance, with or without a buffer zone. 

In addition, curve fits of the analytical data yield equations which express the heat flows 
and degraded radius as functions of the input variables. 



INTRODUCTION 


High-performance multilayer insulation has been proposed to enable long-term 
storage of liquid hydrogen in space. The extremely low heat flux of such an insulation 
system is realized when the heat transfer normal to the shields is primarily by radiation. 
Such a system can be approximated by a large number of highly reflective shields (e. g. , 
aluminized polyester film) separated by low thermal conductivity spacers (e. g. , silk 
netting) . 

The thermal performance of multilayer insulation, however, can be seriously de- 
graded in the area where a penetration such as a plumbing line or tank strut passes 
through the insulation blanket. To prevent any direct radiation from impinging on the 
tank itself, the individual shields must butt up against the penetration. A cross-sectional 
view of a five-shield multilayer insulation system and penetration is shown in figure 1(a). 
Using the minimum thickness of aluminum considered adequate to ensure a low emissiv- 
ity surface, together with a low thermal conductivity substrate (polyester film), reduces 
the lateral thermal conductivity by a factor of 100 over the shield of pure aluminum. 
However, even the use of a minimum thickness of aluminum with its high thermal conduc- 
tivity results in a relatively high thermal conductivity for the shield. 

A technique proposed by several investigators (e. g. , ref. 1) for minimizing the heat 
conducted laterally along the shields uses a thermal buffer zone between the penetration 
and the shields. A cross-sectional view of this configuration is shown in figure 1(b). Un- 
like the radiation shields, the buffer zone has the same thermal conductivity in both the 
normal and lateral directions. This conductivity would greatly exceed the effective ther- 
mal conductivity normal to the shields, but be less than the effective thermal conductivity 
of the shields in the lateral direction. Increasing the buffer zone width would decrease the 
lateral heat flow into the penetration. However, this might be offset by an increase in the 
normal heat flow into the tank. 

The primary objectives of the analysis are to determine representative values for the 
direct lateral heat flow to a penetration, the net increase in heat flow into a tank, and the 
radius of the resulting thermally degraded area. The results of the program are intended 
to emphasize the potential magnitude of the problem, and as such can only be used as a 
guide to preliminary design considerations. For instance, a comparison of the no buffer 
zone configurations with those using a buffer zone establishes a range of values for the 
maximum thermal conductivity that a buffer zone material may possess and still warrant 
its use. 

The analytical results obtained represent maximum values since thermal contact re- 
sistance in the form of a finite temperature differential at any junction could not be readily 
incorporated into the program. However, specifying artificial temperature profiles at the 
shield-penetration interface, or combining the results of the buffer region for one partic- 
ular run (set of input data) with the results for the shield region of another run gave re- 
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suits that simulated thermal contact resistance. The correction factors corresponding to 
various values of thermal contact resistance make the analytical results more descrip- 
tive of an actual application. 


ANALYSIS 

Shield- Penetration Model 

The analytical model for a typical 30° angular segment of the insulation system having 
the shields butted against a penetration is shown in figure 2(a) . A grid network of points 
is set up throughout the segment. The temperature profile at the penetration (T 1 , 

r > * 

T , T J is specified, and the temperature of the tank wall is set at 37° R 

(20. 6 K). A temperature at each grid point is calculated by performing a heat balance 
about that point. The analysis considers only conduction in the radial direction and only 
radiation in the normal direction. The directions and modes of heat flows for a typical 
point (T r N _g) are indicated in figure 2(a) . A converged temperature profile is obtained 
only after the temperatures have been adjusted so that a heat balance at each and every 
point is obtained. The heat flows conducted along each shield into the penetration and the 
radiation exchange between the shield adjacent to the tank and the tank itself are then 
readily determined. 

The spacing between consecutive grid points in the radial direction could take on 
three different values for any particular run. The smallest grid increment (usually 
0. 5 in. or 1. 3 cm) was used in the vicinity of the penetration where changes in the shield 
temperatures were most abrupt. The intermediate grid increment (1.0 in. or 2. 5 cm) 
followed by the largest grid increment (2.0 in. or 5.1 cm) were conveniently set as mul- 
tiples of the smallest grid increment. Several trial runs for 10 shields indicated no sig- 
nificant change in the temperature profiles, and thus the resulting heat -transfer rates, 
when the smallest grid was reduced to 0. 25 or 0. 125 inch (0. 65 or 0. 33 cm) . 

An acceptable solution for any given number of shields occurred when the calculated 
temperature profile at the edge of the model approached the undisturbed temperature pro- 
file (T und j, T und 2 > • • •> T und ‘ The undisturbed temperature profile is that profile 
occurring in a section of basic insulation without the penetration (i.e. , no lateral heat 
transfer) . The undisturbed temperature profile at the edge ensured that no lateral heat 
transfer occurred at this point. Hence, the entire thermally degraded area was located 
within the model. 

As stated previously, the analysis considers heat transfer normal to the shields by 
radiation only, and laterally by conduction along the shields. For a given shield, the 
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resistance, TCR, 
percent 



tb) Lateral heat flow to penetration. 



MMane of penetration Radius 

(c) Normal heat flux to tank wall (Q/A) norm . 
Figure 2. - Shield-penetration model. 
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conduction between two adjacent points is the familiar expression 


Q cond “ ^cs, S ~ 

(All symbols are defined in appendix A.) For radiation, the particular temperature at a 
given grid point is considered to extend over the entire incremental radiation area 
( A ^rad^ surrounding this point: 

^rad “ ^rad^i, i-1 (^i, r " ^i-1, r) 

Thus, as far as radiation is concerned, the temperature profile on any shield is a series 
of step changes in temperature. 

The primary assumption in this analysis was that the spacer contribution to conduc- 
tion in either direction is negligible. This assumption is partially validated by the ex- 
perimental results (ref. 2) which show that, for an ideal installation using silk netting as 
the spacer material, conduction contributes only 20 percent of the overall heat transfer 
in the normal direction. This increase over pure radiation heat transfer can be simu- 
lated approximately by increasing the shield emissivity. 

The primary reason for the assumption is, however, the lack of experimental data to 
enable a valid separation of conduction and radiation on a shield to shield basis. Conduc- 
tion may not be a constant value between consecutive shields, but may be highly dependent 
on the temperature profile pure radiation would tend to establish. A theoretical calcula- 
tion of the conduction term would require knowledge of the thermal conductivity against 
temperature and effective cross-sectional areas of the silk netting. These terms would 
be even more difficult to define for lateral conduction into the penetration. The final ob- 
stacle would involve determining how thermal contact resistance might decrease the con- 
duction term. This thermal contact resistance would occur where the shields touch the 
spacer as well as where the spacer could butt against the penetration. 

The real question is what difference the mode of heat transfer in the normal direction 
has on the results of the analysis. As long as the shield emissivity can be changed to 
vary the overall normal heat flow to simulate some conduction, the results appear valid. 
Using pure radiation heat transfer, as opposed to a combination of radiation and conduc- 
tion (if a valid relation was available), may only tend to establish somewhat higher tem- 
peratures on the shields adjacent to the tank. The temperatures on the shields farthest 
from the tank will undoubtedly remain approximately the same for each case. Since these 
shields transfer the majority of the lateral heat, the total lateral heat conducted into this 
penetration remains essentially the same for each case. This is, of course, the param- 
eter of primary concern. 
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The analysis makes two additional assumptions. First, the radiation view factor be- 
tween corresponding incremental areas of consecutive shields is one. This is valid be- 
cause of the relatively high length to gap ratio characteristic of infinitely long parallel 
plates. A shield spacing of 0.020 inch (0.051 cm) and a grid size of 0. 5 inch (1.3 cm) 
offers a minimum length to gap ratio of 25. The second assumption assumes that shield 
emissivity is not a function of temperature. This assumption is made because experi- 
mental data for variations in emissivity of aluminized polyester film with temperature 
were not available. Because of the relatively high temperatures even in the shield adja- 
cent to the tank, it can be shown that possible variations in emissivity will not signifi- 
cantly change the normal heat flux. 

The direct lateral heat flow to the penetration (Qp j a ^_) is defined as the sum of the 
individual shield heat flows (Qp .) : 


N 


^p , lat X/ ^P , i 
r=l 


The individual heat flows should increase for those shields having a higher undis- 
turbed temperature, as shown in figure 2(b) . 

The net increase in heat flow into the tank (AQp) must consider the decreased normal 
radiation in the vicinity of the penetration. The anticipated normal heat flux resulting 
from the reduced temperatures of the shield adjacent to the tank is shown in figure 2(c) . 
The value of AQp was defined as the total quantity of heat transferred through an arbi- 
trarily defined area surrounding the penetration minus that quantity transferred 

through the same area of undisturbed insulation: 


iQ p- Q p,lat + ^ 


Q, 


norm 


dist 


norm, und 


dist 


where is defined by the radius where the normal heat flux approaches 99 percent of 

the undisturbed heat flux (Q/A) norm und (see fig. 2(c)). 

The analytical program assumes a perfect thermal short between the shield edges 
and the penetration (i.e. , no thermal contact resistance (TCR) is considered). This 
means that whatever temperature profile is imposed on the penetration is simultaneously 
imposed on the shield edges butting against the penetration. If the penetration is at liquid 
hydrogen temperature (T p j = T p 2 ~ ' • • T p N = 37 ° R ( 20 - 6 K))> Qp is a maximum 
for which TCR is defined as zero. If the temperature profile imposed on the shield 
edges corresponded to the individual undisturbed shield temperatures (Tp . = T. und , 
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Heat flow 


Penetration-^ 


(a) Shield-buffer-penetration grid. 


Ul IU, (. 

T und, 1 -^ 


Thermal contact re- 
sistance, TCR, percent 



(b) Normal heat flux to tank wall (Q/A) norm , 
Figure 3. - Shield-buffer-penetration model. 










where i = 1, 2, . . . , N), Qp lat would be zero - a condition for which TCR is defined 
as 100 percent. 

The temperature difference across this interface, corresponding to a given value of 
TCR, determines the resulting heat flow across this interface. The transfer of higher 
lateral heat flows, associated with those shields exposed to higher undisturbed tempera- 
tures, requires a proportionately higher AT to be imposed at this interface. The analy- 
sis uses a given value of TCR to reduce the maximum available AT for each shield by 
a fixed percentage. All shield edge temperatures (Tp .) for a given value of TCR can 
then be calculated in terms of their respective undisturbed temperatures (T un( j .) and 
liquid hydrogen temperature (T^) . The only simplification here is that TCR remains 
constant for each shield of a given penetration insulation configuration.* 

T _ t 

TCR(percent) = 100 — Bi* 

^und, i ~ ^lh 

Shield - Buffer Zone - Penetration Model 

The analytical model for a typical 30° angular segment of the insulation system using 
a thermal buffer zone between the shields and the penetration is shown in figure 3(a) . A 
grid network of points is set up in the buffer zone in addition to the previously described 
grid in the shield region. The buffer zone points are located in the same planes as the 
shields, and they generally had a radial spacing equal to the shield spacing. Again the 
temperatures are adjusted until a heat balance at each and every point is obtained. The 
heat flows out each surface of the insulation segment are then readily determined. 

The emissivity of the top surface of the buffer zone was set equal to the emissivity of 
the shield in order to compute the radiation received by this surface. This can be dupli- 
cated in an actual application by attaching a shield section to the top of the buffer zone. 
This radiation and the heat conducted laterally from the shields are now redistributed into 
two components - that conducted laterally into the penetration (Qg j a ^) and that conducted 
into the tank in the normal direction (Qg norm ) • The total heat transferred out of the 
buffer zone (Qg can then be expressed as 

^B, tot “ ^B, lat + ^B, norm 

The net increase in heat flow into the tank (AQp g) must take into consideration the 
anticipated reduced normal heat flows in the shield region adjacent to the buffer zone as 
shown in figure 3(b). When an arbitrarily defined thermally degraded area (A^.^ g) is 
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used, AQ^g 


can be expressed as 


AQ p, B " Q B, tot + 



dist 


where A digt is again defined by the radius where the normal heat flux approaches 
99 percent of the undisturbed heat flux. 

The analysis assumed a perfect thermal short existed at all interfaces. An actual 
application would have some thermal contact resistance (TCR) at each of the various in- 
terfaces. For convenience, this analysis considered all of the TCR to occur at the buf- 
fer zone - shield interface. The error incurred in not considering TCR at the buffer 
zone - penetration or buffer zone - tank interfaces is only a matter of definition. The ef- 
fect of TCR is to reduce heat flows through the buffer zone by imposing finite tempera- 
ture differences at the interfaces considered. Whether or not the same decrease in heat 
flow is obtained by a series of temperature drops or by one temperature drop only changes 
the relation of TCR against AT at the interface in question. 

Any specific thermal conductivity of the buffer zone material (kg ay g) established a 
specific temperature profile at the buffer zone - shield junction. Various values of 
kr, were used with the same geometric configuration to generate different tempera- 
ture profiles, and their corresponding lateral heat flows across this interface. Since 
this heat flow must be transferred through the buffer zone, a single interface temperature 
(T. t representative of a given profile was used to express the relation between this 
profile and its associated lateral heat flow: 


^S, lat 00 k B, avg^int, avg ” 


The lateral heat flow obtained with a low 
representative of the shield edge profile. 


k B, avg buffer zone y ields a certain T int, avg, S 
This same lateral heat flow can be used with a 


higher kg ayg buffer zone to give a different T int ay g g representative of the buffer 
edge temperature profile: 


Ss,lat* 



int, avg, S 



int, avg, B 



The theoretical temperature difference (T- nt ay g g - T- nt ayg B ) divided by the maxi- 
mum potential temperature difference available (T und g - T int ayg g) defines the value 
of TCR. A relation between TCR and the decreased heat flows is then readily deter- 
mined. 
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Equations and Boundary Conditions 


The equation governing the two-dimensional heat transfer in the buffer region (see 
fig. 4) is 


— [k B (T) — 1 + i-[k B (T) — 1=0 

9r L 9r J 9z L 9z -1 


( 1 ) 


and in the shield region the equation for the top shield is 


W ^P ks ^dr"] +ra [^ rs > NT rs " ^ F rs,N + F N,N-l^ T N,r + F N,N-l T N-l,rJ 0 


and for all other shields i, it is 


> ( 2 ) 


W 


+ ~ ( F i,i + 1 + F i,i-l> T f,r- F 


i,i-l T i-l,rJ 


= 0 


J 


These are the two differential equations to be solved. The solutions for both the no- 
buffer zone and buffer zone models must satisfy the following boundary conditions. Along 
the pipe wall at r = R (fig. 4) the temperature is prescribed: 

r 


T < R p, z> - T p 


(3) 


where T is the local pipe temperature. 
P 


Penetration 

t 


-Ri, 


Pene- 
tration^ 
wall — ^ 


int 


Buffer region 


^dist 


-Interface 


Shield region 


-Tank wall 


Figure 4. - Cross-sectional view of shield-buffer-penetration model. 
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The temperature along the tank surface at z = 0 is constant, so 


T(r, 0) = T m 


( 4 ) 


Since there is no heat flow across the line r = R^.^, 


On the surface of the top shield at 
shield region 


z 



z max 


= 0 (5) 

dist 

there is radiation from space, so that in the 


Q 

^z=z 


= F 


max 


rs 


N"( T rs- T N,r) 


( 6 ) 


Boundary condition equations (7) and (8) apply to the buffer zone model only. 


ST 

Sz 


= F_ 


! ’rs,N CT ( ri 


rs 


2=2 


max 



( 7 ) 


where T _ is the source temperature. There must be continuity of heat flow at the in- 
r s 

terface (r = between the buffer and shield regions; hence, 



( 8 ) 


The equations were solved numerically on a digital computer using a finite-difference 
line-by-line overrelaxation method. Details of the numerical method of solution are given 
in appendix B . 


Program Variables 

The radiation shield simulated in the analysis was composed of a 1/4-mil- 
(0. 00064-cm-) thick polyester film substrate having 600 A of aluminum vapor deposited 
on each side. This was the only aluminum coating thickness used, since experimental 
data (ref. 2) indicated this as a minimum thickness to assure a low emissivity. 
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40 120 200 280 360 440 520 

Temperature, °R 



40 80 120 160 200 240 280 


Temperature, K 

Figure 5. - Comparison of theoretical and experimen- 
tal values of lateral thermal conductivity of a single 
radiation shield as function of temperature. 


A shield emissivity of 0.024 was used for the majority of the data runs. However, 
values of emissivity up to 0.096 were used to determine what effect changing the normal 
radiation had on the magnitude of heat lost to the penetration. Increasing the normal ra- 
diation would be one way of simulating some conduction heat transfer normal to the 
shields that otherwise could not be considered. However, the lateral heat flows down the 
shields may be slightly biased, since changing the emissivity does not alter the resulting 
temperature profiles on the shields. 

The analysis considered two distinct values (experimental and theoretical) for the 
thermal conductivity of the shield. As indicated in figure 5, a significant discrepancy be- 
tween the two occurs over the lower temperature range. The theoretical value is the cal- 
culated effective weighted value of 0.00025 inch (0.00064 cm) of polyester film combined 
with the 1200 A of aluminum, using the bulk thermal conductivity of 99.9-percent pure 
aluminum. The experimental value as determined by Lockheed (ref. 3) indicates the pos- 
sibility that the thermal conductivity of thin films at low temperatures may deviate sig- 
nificantly from the bulk value. The difference is attributed to the effect that the electron 



13 



mean free path has on the thermal conduction properties. Theoretically this effect be- 
comes extremely significant as the electron mean free path approaches the film thick- 
ness. 

The number of shields used was varied from 5 to 160 for the shield penetration 
model, and from 5 to 80 for the thermal buffer zone configurations. Increasing the num- 
ber of grid points in the analytical model for the higher number of shields increased the 
computer program time to the point where the cost in running more than 80 shields with a 
buffer zone became exorbitant. Since only radiation occurred in the normal direction for 
the shield penetration model, the individual shields could be clustered in groups of five. 
The emissivity of each group now acting as a single shield was adjusted so that the normal 
radiation in the undisturbed section of insulation remained unchanged. The effective 
thickness of the clustered shield was increased by a factor of five to determine the lateral 
conduction for the total number of shields simulated. 

Clustering the shields was done for all shield-penetration configurations using 40 or 
more shields. Sufficient points were run comparing clustered to nonclustered configura- 
tions to establish the validity of the technique. 

The shield spacing used in the shield-penetration model was immaterial, since heat 
transfer in the normal direction was by radiation only. When using the buffer zone con- 
figuration, however, the shield spacing was set at 0.020 and 0.010 inch (0.051 and 
0.025 cm). This simulated the "best" and "worst" cases, respectively, for an actual 
application. 

The penetration diameter was generally set at 2 inches (5. 1 cm). However, suffi- 
cient runs were made to determine any performance changes due to changes in the diam- 
eter, both with or without the buffer zone. 

The radiating source temperature was generally set at 520° R (289 K) to simulate the 
extreme case and to maximize the various heat flows being compared. Sufficient runs 
over the range of 200° to 520° R (111 to 289 K) were made to establish performance 
curves. The emissivity of the radiating source was set at 0. 8 and the emissivity of the 
tank surface was set equal to the shield emissivity. 

The average thermal conductivity of the buffer zone material was varied from 0.0004 
to 0. 0095 Btu per foot per hour per °R (0. 00069 to 0. 0164 J/(m)(sec)(K)) . The higher 
value is representative of a 2-pound-per-cubic-foot (32.04-kg/m ) density polyurethane 
foam. The lower value typical of glass wool fibers represents the lowest practical value 
attainable . 

In all buffer zone configurations the buffer zone height was set equal to the insulation 
thickness. The buffer zone width was varied from near zero to twice the buffer zone 
height. 
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RESULTS AND DISCUSSION 


Case 1 - Shields Thermally Shorted to Penetration 

The temperature profiles on the 1st, 10th, and 20th shields of a 20-shield system 
thermally shorted to a 2-inch- (5. 1-cm-) diameter penetration at liquid hydrogen tem- 
perature are shown in figure 6. As shown the shield temperatures are severely degraded 



0 4 8 12 16 20 24 28 32 

Distance, in. 



Distance, cm 


Figure 6. - Temperature profiles for 1st, 10th, and 20th shields of a 
20-shield system thermally shorted to a 2-inch- (5. 1-cm-) diameter 
penetration as function of distance from penetration. Radiating source 
temperature, T rs , 520° R (289 K); shield thermal conductivity, k$, ex- 
perimental; shield emissivity, € s , 0.024. 


15 



in the vicinity of the penetration. The severity of the temperature profile (AT/Ar ratio) 
determined the magnitude of heat conduction along any particular shield into the penetra- 
tion. This ratio increases for shields farther away from the tank and reaches a maximum 
for the top shield. 

The increasing ratio for shields nearing the top results from two independent effects. 
First the shields nearing the top are exposed to increasingly higher undisturbed shield 
temperatures, thus increasing the potential AT available. The undisturbed shield tem- 
perature is that temperature which exists on the shield when the lateral heat conduction 
approaches zero. Secondly, the distance required for a given shield to recover to its un- 
disturbed temperature increases for shields nearer the tank wall. This has the effect of 
increasing the effective heat path. This trend is independent of the number of shields 
comprising the insulation system as shown in figure 7. 

A qualitative explanation of this trend considers a heat balance about the top shield. 
The total quantity of heat transferred through an area of insulation surrounding a penetra- 
tion is dependent solely on the temperature profile of the top shield. This profile deter- 



Figure 7. - Radius about a 2-inch - (5. 1-cnH 
diameter penetration for the top and bottom 
shields to reach a specific temperature dif- 
ferential from the undisturbed temperature 
for that shield as function of number of 
shields in the system. Radiating source 
temperature, T^, 520° R (289 K); shield 
thermal conductivity, k<r, experimental; 
shield emissivity, € s , 0.024. 
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mines the net exchange of radiation between the radiating source and the top shield, as 
expressed in 


Siorm, rad AAcrF rs, rs " T N, r) 


V 


If the conduction term along the top shield were zero, all radiation received by the top 
shield would necessarily be passed on to the shield below. Considering F rg ^ equal to 
F N,N-1 results in the temperature relation 





•) 


Table I illustrates how small changes in T N.r cause significantly larger changes in 
T n _i r for a constant radiating source temperature of 520° R (289 K). Thus for tem- 
peratures on the top shield less than 437° R (243 K) , the shield below could not attain 
the required temperature to keep the net exchange approximately the same. To compen- 
sate for this, the incremental areas on the second shield must be larger than their cor- 
responding areas on the first shield. This same argument can be used to show why the 
degraded area for any other shield must be less than the degraded area of the shield be- 
low it and greater than the degraded area of the shield above it. 


TABLE I. - SHIELD TEMPERATURE 
COMPARISON FOR PURE RADIATION 
HEAT TRANSFER 


[Radiating source temperature, 520° R 
(289 K) . ] 


Temperature 

T N, r 

t n 

-l,r 

°R 

! 

K 

°R 

K 

500 

278 

477 

265 

480 

268 

426 

237 

460 

256 

358 

199 

440 

244 

207 

115 

437 

243 

69 

38 
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0 30 60 90 120 150 180 210 240 270 

Distance from penetration, cm 

Figure 8. - Normal heat flux as function of distances from a 2-inch- (5. 1-cm-) diameter penetra- 
tion for various number of shields. Radiating source temperature, T. s , 520° R (289 K); shield 
thermal conductivity, k s , experimental; shield emissivity, €$, 0.024. 



Figure 9. - Radius of thermally degraded area (Q norm = 0. 99 
Q n0 rm und* arouncl 2-inch- (5. 1-cm-) diameter penetration as 
function of number of shields. Radiating source temperature, 
T rs , 520° R (289 K). 
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Heat flow, Btu/hr 




(b) Net increase in heat flow, AQ p . 

Figure 10. - Heat flow into tank fora 2-inch- (5. 1-cm-) diameter 
penetration as function of number of shields. Radiating source 
temperature, T rs> 520° R (289 K). 
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Number of shields 

(b) Shield thermal conductivity, k$, experimental. 

Figure 11. - Net increase in heat flow into tank for various diam- 
eter penetrations as function of number of shields. Radiating 
source temperature, T rs , 520° R (289 K). 
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Temperature, K 


Figure 12. - Net increase in heat flow into tank for 2-inch- (5. 1-cm-) 
diameter penetration as function of undisturbed top shield tempera- 
ture. Shield thermal conductivity, k s , experimental; shield emis- 
sivity, € 5 , 0.024. 


The degraded temperatures on the shield adjacent to the tank decrease the normal 
radiation from this shield to the tank as shown in figure 8 . The end point on each curve 
represents the arbitrarily defined disturbed radius where the normal radiation component 
equals 0.99 of the undisturbed heat flux. The disturbed radius sets the degraded area 
used for determining the net heat flow into the tank. 

The radii defining the thermally degraded area for various number of shields for two 
different shield thermal conductivities and for various shield emissivities are shown in 
figure 9. In general, doubling the number of shields tends to double the radius of the 
thermally degraded area for any set of variables. As the shield emissivity increases, 
the radius decreases. Changing the normal component by radiation does not significantly 
change the heat conducted laterally along the shield. Thus, any particular shield does 
not have to draw on as large an area to essentially transfer the same amount of heat to 
the penetration. 

The direct lateral heat flow to the penetration Qp ^ and the net increase in heat 
flow into the tank AO as a function of the number of shields and shield emissivities are 
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shown in figures 10(a) and (b) , respectively. In general, doubling the number of shields 
increases the respective heat flows by approximately 50 percent. Doubling the number of 
shields doubles the effective number of individual heat paths. However, this potential in- 
crease is partially offset by the fact that doubling the number of shields also increases 
the effective heat path for any given shield (i.e. , the disturbed radius increases). For 
any given set of conditions AQ^ is approximately equal to 0. 67 Qp lat . 

The effect of various diameter penetrations on AQp as a function of the number of 
shields for both the theoretical and experimental values of shield thermal conductivity is 
shown in figure 11. Increasing the penetration diameter has less effect on AQ p on a 
percentage basis as the number of shields is increased (i.e. , the curves converge). 

This trend indicates that for the low number of shields the shield cross-sectional area 
is more influential than the number of individual heat paths in determining the lateral heat 
flows. For the high number of shields the reverse is true. 

The effect of changing the undisturbed top shield temperature on the net heat flow into 
the tank (AQp) is shown in figure 12. Various top shield temperatures are obtained by 
exposing the shield to different radiating source temperatures . Since radiation heat 
transfer is proportional to the fourth power of the temperature, the curves in figure 12 
are not linear. The data from figure 12 are normalized in figure 13 to a reference tem- 
perature of 520° R (289 K) . The normalized curve is independent of the number of 
shields being used. 



40 80 120 160 200 
Temperature, K 


240 280 


Figure 13. - Ratio of net increase in heat flow into tank for a given 
undisturbed top shield temperature to that heat flow when this 
shield temperature is 520° R (289 K). Shield thermal conductivity, 
k s , experimental; shield emissivity, e s> 0.024. 
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Case 2 - Thermal Contact Resistance for Shield Penetration Model 


All of the data presented in Case 1 assumed a perfect thermal short between each 
shield edge and the penetration. An actual application would offer some finite thermal 
contact resistance (TCR) even if the installation were perfect (i.e. , no gap between the 
shield edge and the penetration) . 

The value of TCR as previously defined in the analysis decreases the lateral heat 

flow along each shield by specifying the shield edge temperature T . : 

P> ^ 


TCR = 100 P)1 “L_ 

T und, i " T lh 

The shield edge temperature profiles for various values of TCR for a 20-shield insula- 
tion system are shown in figure 14. Curves of this type can be used to determine a range 



Figure 14. - Shield edge temperature profiles for a 20-shield system as 
function of thermal contact resistance. 
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Figure 15. - Ratio of direct lateral heat flow to a penetra- 
tion and/or net increase in heat flow into the tank to its 
respective thermally shorted value as function of thermal 
contact resistance. 


of values of TCR to use in terms of an anticipated temperature difference between a 
given shield and the penetration. Specific values of thermal contact resistance to use 
will have to await experimental testing. Quite possibly the installation technique and care 
in handling will greatly influence the value of TCR to be determined. 

The ratios of Qp ^ and/or AQp to their respective thermally shorted values as a 
function of TCR are shown in figure 15. The ratios are independent of the number of 
shields, but highly dependent on the thermal conductivity of the shield. This is due to the 
fact that the experimental value of thermal conductivity decreases with lower tempera- 
tures, while the theoretical value does the opposite. 


Case 3 - Use of Thermal Buffer Zone Between Shields and Penetration 

The data presented in this section consider the shields thermally shorted to the buffer 
zone. The buffer zone in turn is thermally shorted to the penetration and tank. All tem- 
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Temperature, 







peratures on the tank surface and penetration are set at 37° R (20.6 K). Unless specified 
otherwise, the following values for the variables are used: 


Shield emissivity, eg 0.024 

Penetration diameter, D^, in. (cm) 2(5.1) 

Shield spacing, AZg, in. (cm) 0.020 (0.051) 

Shield thermal conductivity, kg experimental 

Radiating source temperature, T _, °R (K) 520 (289) 

AO 


The temperature profiles through the buffer zone corresponding to the planes of the 
1st, 6th, 11th, 16th, and 20th shields of a 20-shield system are shown in figure 16. The 
two average thermal conductivities kg av ^ of the buffer zone material (0.0095 and 
0.0004 Btu/(ft) (hr) (°R) ; 0.0164 and 0.00069 J/(m)(sec)(K)) represent the extremes for 
appropriate insulation materials . In each case the buffer zone width was arbitrarily set 
equal to the buffer zone height of 0.40 inch (1. 02 cm). 

As would be expected, the lower k-o yields the higher temperature profile at the 

r>, avg 

shield - buffer zone interface. As the temperature profile at tins interface increases, the 
quantity of heat transferred laterally from the shields to the buffer zone must decrease. 
The optimum profile at this interface would be the undisturbed temperatures of the indi- 
vidual shields. 

The effect of the buffer zone width and thermal conductivity on the heat transferred 
through the buffer zone for a 20-shield system is shown in figure 17. The two buffer zone 
heights of 0. 20 and 0. 40 inch (0. 51 and 1.02 cm) represent shield spacings of 0. 010 and 
0.020 inch (0. 025 and 0.051 cm), respectively. In each configuration, increasing the 
buffer width decreased the lateral heat flow into the penetration (Qg j at ), but increased 
the normal heat flow into the tank (Qg norm ) • The optimum geometric configuration is 
that configuration which minimizes the total heat flow into the tank (Qg ^ + Qg norm ) • 

This minimum heat flow occurred at a buffer width approximately equal to the buffer 
height, and at a width three -fourths of the height for kg ay g = 0. 0004 Btu per foot per 
hour per °R (0.00069 J/(m)(sec)(K)) . For kg ayg = 0.0095 Btu per foot per hour per °R 
(0. 0164 J/(m)(sec)(K)), this minimum heat flow occurred at a buffer width varying from 
one -half to three -fourths of the buffer height. However, the net increase in heat flow into 
the tank for a buffer material having k„ = 0.0095 Btu per foot per hour per °R 

.Dj avg 

(0.0164 J/(m)(sec)(K)) was approximately that obtained without a buffer zone. (This is 
discussed in more detail in conjunction with fig. 20.) All new data presented will approxi- 
mate optimum values for acceptable buffer zone materials by using a square configura- 
tion. 

Using a shield spacing of 0.010 or 0.020 inch (0.025 and 0.051 cm) (B w = B^ = 0.2 
and 0.4 in. (0.51 and 1.02 cm), respectively, for a 20-shield system) has very little ef- 
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(a) Buffer zone height, B h , 0.2 inch 
(0. 51 cm); average buffer zone 
thermal conductivity, kn a 
0. 0095 Btu per foot per nbur per °R 
(0. 0164 J/(m)(sec)(K)). 



(b) Buffer zone height, B h , 0.4 inch 
(1.02 cm); average buffer zone 
thermal conductivity, kg avg , 
0.0095 Btu per foot per nbur per °R 
(0.0164 J/(m)(sec)(K)>. 



0 .2 .4 .6 .8 1.0 

Width, cm 

(c) Buffer zone height, B h , 0.2 inch 
(0.51 cm); average buffer zone ther- 
mal conductivity, kg avg , 0.0004 
Btu per foot per hour per °R 
(0. 00069 J/(m)(sec)(K». 


0 .4 .8 1.2 1.6 2.0 


(d) Buffer zone height, B h , 0.4 inch 
(1.02 cm); average buffer zone ther- 
mal conductivity, kg avg , 0.0004 
Btu per foot per hour per °R 
(0. 00069 J/(m)(sec)(K)). 


Figure 17. - Effect of buffer zone width on heat flows through a buffer zone ther- 
mally shorted to a 2-inch - (5. 08 -cm-) diameter penetration and a 20-shield 
system. Radiating source temperature, T rs , 520° R (289 K); shield emissivity, 
0.024; shield thermal conductivity, k$, experimental. 
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Figure 18. - Effect of shield spacing on net increase in 
heat flow into tank as function of number of shields. 
Radiating source temperature, Tj- S , 520° R (289 K); 
shield emissivity, 0.024; shield thermal conduc- 
tivity, k s , theoretical; penetration diameter, Dp, 

2 inches (5. 1 cm). v 
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Figure 19. - Effect of shield thermal conductivity on net 
increase in heat flow into the tank as function of num- 
ber of shields. Radiating source temperature, T rs> 

520° R (289 K); shield emissivity, c s , 0.024; shield 
spacing, Az, 0.020 inch (0.051 cm); penetration diam- 
eter, Dp, 2 inches (5. 1 cm). 
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Figure 21. - Net increase in heat flow into the tank as function of average thermal conduc- 
tivity of buffer zone. Number of shields, 20; radiating source temperature, T rs , 520° R 
(289 K); shield spacing, Az, 0.020 inch (0.051cm); penetration diameter, D p , 2 inches 
(5.1cm). 
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Figure 22. - Ratio of net increase in heat flow into the tank for a 
given undisturbed top shield temperature to heat flow when this 
shield temperature is 520° R (289 K). Number of shields, 20; 
shield thermal conductivity, k$, experimental; shield emissivity, 
0.024; average buffer zone thermal conductivity, kn avq , 
0.0095 Btu per foot per hour per °R (0. 0164 J/(m)(sec)(K)l. 


Penetration Average buffer zone thermal 
diameter, D 0 , conductivity, kg avg , 
in. (cm) Btu/(ft)(hr)(°R) (Jj(m)fsec)(K)) 
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Figure 23. - Net increase in heat flow into the tank using a 
buffer zone around various diameter penetrations as function 
of number of shields. Radiating source temperature, T fS , 
520° R (289 K); shield thermal conductivity, k$, experimental 
shield emissivity, 0.024. 


feet on AQp, g. This is shown in figure 18 for kg aV g = 0. 0095 and 0.0004 Btu per foot 
per hour per °R (0. 0164 and 0.00069 J/(m)(sec)(K)) and the theoretical value of kg. The 
significant change in AQ^ occurs with changes in kg ay g for a given number of shields. 
The effect of shield thermal conductivity on AQp g becomes less significant as 


k R decreases. As shown in figure 19, the use of either experimental or theoretical 
o, avg 

results in the same AQp g when kg ayg = 0.0004 Btu per foot per hour per R (0.00069 
J/(m)(sec) (K)) . 

Doubling the number of shields increases the AQp g by approximately 50 percent 
for both values of kg ay g as shown in figure 20. This relation is exactly the same as 
was previously shown for the shields butted directly to the penetration. The comparable 
data from figure 10(b) for the shield-penetration configuration are replotted in figure 20. 
This curve indicates that use of a buffer zone for k R = 0.0095 Btu per foot per hour 
per R (0.0164 J/(m)(sec)(K)) is thermally inferior to the no-buffer zone configuration. 
For buffer zone thermal conductivities in this range, the decreased lateral heat flow was 
offset by the increased normal heat flow through the buffer zone. Since all known conven- 
tional buffer zone materials are isotropic, the analysis considered the thermal conduc- 
tivity equal in both the lateral and normal directions. Thus, a reduction of the lateral 
heat flow could be accomplished only by an increase in the normal heat flow. As shown, 
there is a limiting value of thermal conductivity that a buffer zone can have to still war- 
rant its use. 

The effect changing kg ay g has on AQp g for various shield emissivities is shown 
in figure 21. The value of eg used becomes insignificant only at low values of kg ay g. 
At the higher values of kg ay g, the increase in eg creates a higher potential quantity of 
heat to be transferred laterally into the buffer zone. 

The ratio of AQp g for a given undisturbed top shield temperature to the AQp g 
when this top shield temperature is 520° R (289 K) is shown in figure 22. Again the curve 
is nonlinear since radiation heat exchange is proportional to the fourth power of the tem- 
perature. 

The effect of changing the penetration diameter on AQp g is shown in figure 23. As 
shown previously for the shield-penetration model, the curves again tend to converge at 
the high number of shields. 


Case 4 - Thermal Contact Resistance for the Buffer Zone 

All of the results presented for Case 3 assume that the buffer zone is thermally 
shorted to the edges of the shields as well as to the penetration and tank surface. This 
results in maximum heat flows through the buffer zone, and thus eventually into the tank. 
To adjust these results for an actual application, a correction factor expressing thermal 
contact resistance has been generated. 
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Figure 24. - Temperature differentia! across buffer zone - shield 
junction in plane of 11th shield of 20-shield system as function of 
thermal contact resistance. 


Thermal contact resistance has been considered at the buffer zone - shield interface 
only. This means that any decreased heat flow through the buffer zone will be realized 
by reducing the lateral heat flow across this interface. This is valid since the data show 
that approximately 90 percent of the total heat flow through the buffer zone does in effect 
cross this interface. The other 10 percent is by direct radiation received by the top sur- 
face of the buffer zone. 

The resulting AT at the buffer -shield interface to simulate various values of 
thermal contact resistance is shown in figure 24. At TCR = 0, a thermal short exists; 
and for TCR = 100 percent, the shield edge temperature is the undisturbed temperature 
for that shield. Experimental measurements of the existing AT across the interface 
will enable calculation of a range of actual TCR values to use in future analyses. 

The ratio of AQp g to its thermally shorted value as a function of TCR is shown in 
figure 25. Various values of k-r, as well as two different number of shields were 
used to generate the data. As indicated, the results show a well defined relation. 
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Figure 25. - Ratio of net increase in heat flow into the tank for a 
given buffer zone configuration to its respective thermally shorted 
value as function of thermal contact resistance. Shield thermal 
conductivity, k s , experimental. 


General Comparisons and Conclusions 


The data presented in the preceding four cases can be used to estimate the degree of 
thermal degradation attributed to several representative penetration configurations. 
However, these results cannot be clearly presented on a single figure because of the many 
variables involved. To avoid scanning numerous figures to estimate the performance, 
curve fits of the data obtained in both Cases 1 and 3 have been made. These curve fits 
are equations which approximate the effect each variable has on the heat flows and de- 
graded radius. For the case of the shields thermally shorted to the penetration the fol- 
lowing are true: 
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For the configuration of the buffer zone between the shields and the penetration the fol- 
lowing equations are true: 
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where D is expressed in inches, T,. N in °R, and k R _ in Btu per foot per hour 

P ID, civ g 

per °R. 

The net increase in heat flow into the tank is approximately two-thirds of the direct 
lateral heat flow to the tank and/or penetration in either configuration. The heat flows 
are extremely temperature dependent and least affected by changes in emissivity. The 
disturbed radius in each configuration increases when either the emissivity or radiating 
source temperature decrease. 

The previous relations used the experimental values of the lateral thermal conduc- 
tivity of the shield against temperature (ref. 3). The shield simulated a 1/4-mil 
(0.00064-cm) polyester film aluminized on each side to a thickness of 600 A; thus, the 
lateral heat conduction path was minimized. The results in graphical form in the body of 
the report use a theoretical value of thermal conductivity as well. The theoretical value 
is obtained by using the bulk properties of aluminum, and quite possibly does not apply to 
thin films. Therefore, the results using the theoretical value may in reality apply to a 
much thicker coating, and thus may only serve to indicate the magnitude of differences 
for various coating thicknesses. 

The maximum thermal conductivity that a buffer zone material may have and still be 
thermally competitive with the no -buffer zone configuration can be determined by equating 
the previous expressions for AQp g and AQp g. Table II lists the wide range of maxi- 
mum values of k R that the buffer zone may have. 

The validity of these comparisons may depend on the experimental value of thermal 
contact resistance for an actual application. Because of this and other factors, as well 
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TABLE n. - BUFFER ZONE THERMAL CONDUCTIVITIES GIVING SAME 


THERMAL PERFORMANCE AS THEIR CORRESPONDING 
NO-BUFFER ZONE CONFIGURATION 


Number 

Shield 

Penetration 

Temperature, 

Average thermal conductivity 

of 

emissivity, 

diameter, 

t s.n 

of buffer zone material, 

shields, 

€ S 

D 

P 



k B,avg 

N 


in. 

cm 

°R 

K 

Btu/ (ft) (hr) (°R) 

j/(m)(sec) (K) 

20 

0.024 

2 

5.1 

520 

289 

0.0076 

0.0131 

80 

. 024 

2 

5.1 

520 

289 

.0054 

.0093 

20 

.072 

2 

5.1 

520 

289 

.0082 

.0142 

80 

.072 

2 

5.1 

520 

289 

.0058 

.0100 

20 

.024 

6 

15.3 

520 

289 

.0055 

.0095 

80 

.024 

6 

15.3 

520 

289 

.0039 

.0067 

20 

.024 

2 

5.1 

400 

222 

.0088 

.0152 

80 

.024 

2 

5.1 

400 

222 

.0063 

.0109 


as different concepts that are not considered in this analysis, the use of a conventional 
thermal buffer zone may not necessarily offer the best overall thermal performance. One 
of these concepts would use an anisotropic buffer zone that has a lower thermal conduc- 
tivity in the normal than in the lateral direction. Thus, the buffer zone width could be 
increased (decreasing the lateral heat flow) without offsetting any apparent gain with a 
high normal heat leak. A potential anisotropic buffer zone would use a separate circular 
segment of basic multilayer insulation immediately around the penetration. Using a sep- 
arate segment interrupts the lateral heat path into the penetration, while still utilizing 
the extremely low thermal conductivity in the normal direction. 

The relation between the heat flows and the equivalent area of basic insulation these 
heat flows represent is shown in figure 26. Thus, the heat flows obtained with the pre- 
vious expressions, either corrected for thermal contact resistance or not, can be related 
to the performance of the overall insulation system. The areas represent maximum val- 
ues since the shield emissivity of 0.024 used in the minimum value that could be antici- 
pated for aluminized polyester film. 

The radius of the thermally degraded area, defined by m = 0. 99 Q norm unc p 

obtained with various values of AQp is shown in figure 27. The aforementioned expres- 
sions for the degraded radius are valid only for the thermally shorted cases . Although 
figure 27 is specifically for the shield penetration model, the results approximate those 
obtained with the buffer zone configurations. 
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Figure 27. - Radius of thermally degraded area as function of net increase in heat flow 
into the tank. Shield emissivity, €§, 0.024. 


CONCLUDING REMARKS 


An analytical investigation was conducted to determine the direct lateral heat flow to 
a penetration the resultant net increase in heat flow into a tank AQ^, and the radius 
R d ist the thermally degraded area created by the presence of a cylindrical pipe or 
tank support penetrating the multilayer insulation on a liquid hydrogen tank. Two differ- 
ent configurations were considered: (1) shield edges butted to the penetration and (2) use 
of a thermal buffer zone between the shields and the penetration. 

The relation between the program variables and the parameters of primary interest 
can be best expressed in the following equations. These equations represent curve fits of 
the analytical data. For the configuration of the shields thermally shorted to the penetra- 
tion the following hold: 
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where 

N number of shields 

€g shield emissivity 

Dp penetration diameter, in. 

Tg ^ undisturbed top shield temperature, °R 

The use of a thermal buffer zone between the shields and the penetration yields the fol- 

lowing: 
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where k R is the average thermal conductivity of the buffer zone material (Btu/(ft) 
(hr)(°R)) . 

The previous equations used the experimentally determined values for the lateral 
thermal conductivity of the shield against temperature (ref. 3) . Equating AQp g (from 
the shield-penetration model) to AQp R (from the shield-buffer -penetration model) shows 
that the buffer zone material must have a thermal conductivity, in general, of less than 
0.008 Btu per foot per hour per °R (0.0138 J/(m)(sec)(K)) to warrant its use. Because of 
thermal contact resistance, as well as the existence or other penetration insulation con- 
cepts, the use of a conventional thermal buffer zone may not necessarily offer the best 
overall thermal performance. An anisotropic buffer zone having a lower thermal conduc- 
tivity in the normal direction than in the lateral direction may be thermally superior. 

The optimum thermal buffer zone configuration (i.e. , minimum net heat flow into the 
tank), in general, has a buffer width approximately three -fourths of the buffer height. 

This was generally true regardless of the number of shields, spacing of the shields, or 
the thermal conductivity of the buffer zone material. 

To make. all of the previous results applicable to a real installation, a correction for 
possible thermal contact resistance at the shield-penetration interface or shield-buffer 
zone interface must be made. A definition of thermal contact resistance for both config- 
urations was established. The ratio of AQp to its thermally shorted value as a function 
of thermal contact resistance is presented. These correction factors can be applied di- 
rectly to the equations for AQp g and AQp R to obtain more representative values of 
heat flows. Unfortunately exact design values for thermal contact resistance will vary 
for each insulation configuration and installation technique. A range of values to use will 
have to be determined by experimental measurement. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, April 17, 1968, 
180-31-08-06-22. 
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APPENDIX A 


SYMBOLS 


A 

f ,2 2 
area, ft ; m 

K 

AA 

2 2 

incremental area, ft ; m 


s/ 

matrix of a system of finite - 



difference equations 


a(x) 

(dN/du/ k_1 ) 

K 

B 

thermal buffer zone 


£f 

submatrices of order N which 



lie along diagonal of 

k 

b(x) 

Nu( k_1 ) - (dN/dU)^ k "^U^ k-1) 



off-diagonal submatrices of 

N 


order N of sf 

NJ 

c i 

coefficients of eq. (B18) defined 



in eq. (B19), 1= 1, 2, . . . , 5 

Q 

D 

diameter, in. ; cm 

AQ 

d 

vector whose components are 



right-hand sides of finite - 

R 


difference equations 

r 

F 

net radiation heat exchange 

Ar 


factor between two gray 


bodies 

T 

G 

00 /ro 

AT 

I 

position number of horizontal 



grid line 

TCR 

J 

position number of vertical grid 

U 


line 


JEND 

total number of vertical grid 

V 


lines 



. . . + T. . 
i,J-i 1,3 


used in 


finite-difference equations at 
interface 


!• j+l + T i, i 


used in 


finite-difference equations at 
interface 

thermal conductivity, Btu/(ft) 
(hr) (°R) ; J/(m)(sec)(K) 

number of shields 

position number of interface 
vertical grid line 

heat flow, Btu/hr; J/sec 

net increase in heat flow into 
tank, Btu/hr; J/sec 

radius, in.; cm 


radial coordinate 


grid spacing in radial direc- 
tion, in. ; cm 

temperature, °R; (K) 

temperature difference, °R; K 

thermal contact resistance 

pseudotemperature (transformed) 
in shield region 

pseudotemperature (transformed) 
in buffer region 
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V 

vector of pseudotemperature 

h 

height 


(transformed) in buffer 

i 

variable position in axial 


region intermediate in the 


direction 


computation 

int 

buffer zone - shield interface 

V 

volume, ft**; m 8 

3 

variable position in radial di- 

w 

shield thickness, in.; cm 


rection 

z 

insulation thickness, in. ; cm 

lat 

lateral 

z 

axial coordinate 

lh 

liquid hydrogen 

Az 

grid spacing in axial direc- 

max 

maximum 


tion, in. ; cm 

N 

top shield position 

Of 

*i, i+l^i+1, r + ^i, i-l^i-1, r 

nom 

nominal 

€ 

emissivity 

norm 

normal direction 

a 

Stefan-Boltzman constant, 

P 

penetration 


0. 17 3X10" 8 Btu/ (ft 2 ) (hr) (°R 4 ) ; 
5.67X10' 8 J/(m 2 )(sec) (K 4 ) 

r 

radial coordinate 


rad 

radiation 

T 

dummy variable 

rs 

radiating source 

A0 

size of angular segment used as 




analytical model, deg 

S 

shield 


block successive overrelaxation 

tot 

total 


parameter 

ts 

thermally shorted 

Subscripts : 

und 

undisturbed 

avg 

average 

w 

width 

B 

thermal buffer zone 

Superscripts: 

cond 

conduction 

(k) 

k*"* 1 iterate 

cs 

cross section 

(m) 

m** 1 iterate 

dist 

disturbed, or thermally de- 

(o) 

initial approximation 


graded 

— 

vector notation 

e 

shield edge 
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APPENDIX B 


NUMERICAL METHOD OF SOLUTION 

The heat flow equations for the buffer (eq. (1)) and shield (eqs. (2)) regions are non- 
linear, and the boundary conditions for these regions are also nonlinear. This precludes 
a closed-form solution so the problem was approached numerically. 

One-to-one transformations were made as follows on equations (1) and (2) to elimi- 
nate the variable conductivity terms. Let 



k B (T)dr 


so that 


and 


^ = wt) " 

dr a dr 


^ = k B (T)^ 

dz a dz 


Substituting equations (B2) and (B3) into equation (1) gives 


(Bl) 


(B2) 


(B3) 


9^V + d^V = 0 


dr‘ 


3z 


(B4) 


Similarly, let 



kg(r)dr 


(B5) 


then 


— = kc(T) — 
dr b dr 


(B6) 
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... J 



Note that in all the equations which follow, the dependence of T on U or V may not be 
explicitly stated. 

Substituting equation (B6) into equations (2) gives 


W— [r— | + ra 


dr\ dr 


[ F rs,N T rs " (F rs,N + F N,N-l )T N,r + F N, N-1 T N-1, r] = 0 

W dr( r d7) + rCT [ Fi > i + lT i+ 1 , r " (F i,i+l + F i,i-l )T i,r + F i, i-l T i-l, rj = 0 
The transformations modify the boundary conditions. At r = R p , 

V(R p>2 )=Vp( T p) 

The values of the dependent variables along the tank surface at z = 0 become 

V = V(T lh ) = 0 


> (B7) 


(B8) 


(B9) 


and 


u = U(T lh ) = 0 


On the adiabatic boundary at r = R^.^, 


dU 

dr 


= 0 


r=R dist 


(BIO) 


(Bll) 


The top surface at z = z max has the same boundary condition (eq. (6)) in the shield re- 
gion as before the transformation 


Q 

^z=z 


max 


' F r S ,N"( T r S - T N,r) 


(B12) 


In the buffer region at z = z max , equation (7) becomes 


(f) = F rs,N"(4- T N, r ) 

' z=z 

max 


(B13) 
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At the interface, the equation for continuity of heat flow (eq. (8)) is now 




(B14) 


The numerical solution was obtained by replacing equations (B4) and (B7) by their 
finite-difference analogues at intersection points of a grid superimposed on the combined 
region of figure 4. These points are shown in figure 28. The solution of the differential 


I-N 



J = NJ 


I I I I 

t t J 1 

I II 


r r~ i 

I T T 


Figure 28. - Typical grid points where solution was obtained. 


I 

I 

I 

T 


J = JEND 


equations satisfies these difference equations except for an error term which vanishes 
when the grid spacing approaches zero. 

As indicated in figure 28, each grid point has an ordered pair of integers (I, J) asso- 
ciated with it. The grid was constructed in such a way that each boundary - including the 
interface (J = NJ) - is a grid line. JEND is the total number of vertical grid lines and N 
is the total number of horizontal grid lines. The horizontals are uniformly spaced as are 
the verticals in the buffer region; the spacing of the verticals in the shield region is not 
necessarily uniform. 

A dual grid was constructed in the buffer region by bisecting the grid lines in that 
region. This resulted in a closed cell about each point of the original grid. The bound- 
aries of the cells are the dashed lines shown in figure 29. 

The finite -difference equations for the pseudotemperatures (V) in the buffer were ob- 
tained as follows: (1) integrating equation (B4) over cells of the dual grid; (2) applying 
the divergence theorem; and (3) replacing the normal derivatives occurring in the bound- 
ary integrals by difference quotients, zero or F rg n ct^T^ s - T^ , depending on the po- 
sition of the cell in question (ref. 4, sec. 6.3). 
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I 


z 


I - N 


I = 0 
J 


<N, 1) 


-interface 


+I+I+I+I4- 

i ± i * 1 ± ! ± i 

tltl+Itlt 
+ I + I+RI + 
tl+ltlflt 

+1 4* J + J-i- J -j- 


J - NJ 



Figure 29. - Dual grid in buffer for numerical solution. 


(N, 1) 



Figure 30. - Detail of grid for numerical solution in buffer. 


The preceding process is illustrated by the following example taken at the point 
(I, J) = (N ? 1) in figure 29. The cell of the dual grid under consideration (shaded in fig. 30) 
has its boundaries (except the upper one) midway between grid lines of the original grid. 
The finite-difference equation at the point (N, 1) is derived by the following steps where 
i = N and j = 1: 
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. Mill Hill || II mu II I III III nil ■Ill am i 


II INI III III I I I II 1 1 ■ II II I II ■■ nil linn 


0 « 


V - v. . 

" hi 


V 


V, 4 


R; 


A0 ** + -iliil hi R. A0 Ar B + -iil +1 




Ar 


B 




Ar 


B 


111 R W l/2) AS 


+ F rs,N"( T J S - T 1 4 ,j) R j Afl Ar B 


Az 

2 


(B15) 


The integration is repeated on all cells that surround the points (I, J) when 1 < I < N 

and 1 < J < NJ - 1; this yields a system of N(NJ - 1) equations in as many unknowns 

V(I,J). Some (NJ - 1) of these equations are nonlinear - namely (N,J), 1 < J < NJ - 1. 

Since this system cannot be solved directly because of the nonlinearity, an iterative 

procedure must be used. Let V^ m ^(I, J) represent the iterate for V(I, J) where 

V^(I,J) is some initial approximation to V(I, J). Equations at the points (N, J), for 

4 

1 < J < NJ - 1, are linearized by replacing the term T (N, J) by the first two terms in 
its Taylor expansion about J) : 



The resulting system of equations linear in can be written in matrix notation 

as 


^y(m) = d 


where 


*-*(. r,; V< m -») 

is a square matrix of order N(NJ - 1) and 

d = d(r,; V^- 1 ) 

is a column vector with N(NJ - 1) components. The block tridiagonal matrix si is given 

t>y 
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I I II I I 


I I I 





a. 


= V^ m '^(N, J) 


1 < J < N J - 1 


are square tridiagonal matrices of order N, and 

*j -*(v v<“>; v<”-«) 


is a diagonal matrix of order N. The iterative method can be described as a line-by-line 
(vertical) overrelaxation scheme, which is carried out as follows. 

An initial approximation to V is made at all points of the buffer region (e.g. , 
V^°\l,J) = V ). Then, V is changed by sweeping through these points an entire vertical 
at one time; that is, the V distribution at grid points on any one vertical is changed 
simultaneously. After a line is changed, the new V vector - denoted by V - is over- 
relaxed: 

~(m) = y(m-l) + w ^(m) _ y(m-l)^ 

where w is the ’’block successive over relaxation parameter” (ref. 4, ch. 4) that con- 
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I-N 


I - N - 1 









\ ' 

1 


— 4 — 


J — 
I 

i 


- 1 — . 


Ar„ 


-Interface 


J = 0 


J - NJ 


Ar s 


Figure 3L - Detail of grid for numerical solution at interface. 


trols the rate of convergence. This scheme is applied to each line (vertical) moving from 
left to right - J = 1 through J = NJ - 1 - before working on the interface (J = NJ) . 

Figure 31 shows a typical cell on the interface surrounding the grid point (I, J) = 

(N - 1,NJ). Since the conductivity functions in the buffer and shield regions are differ- 
ent, the finite -difference equations for the interface were written in terms of the temper- 
ature T . 

Following a procedure analogous to that outlined on page 45 results in the following 
finite -difference equation for the point (N - 1,NJ) where i = N - 1 and j = NJ: 




Az 


Ar 


B 


Vl Ar B - 


1 \ A(,Ar B 


2Az 


+ K i-(l/2),j( T i-l,j - T i,j ) ( r j "I^B, 2Az 


-lAr \ AeAl B 


+ K i,j + (l/2)^ T i,j + l - T 




A0W 

Ar 0 


.4 4 I/ 1 \ Ad Ar S 


+ - T ?,i)( r j + i Ar s)-^r^ 


(B16) 
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where 


-\ 


K i, j-(l/2) ~ k B (T i, Ml/2))- k B 


T i>3-1 + T i>3 


K i,j+(l/2) - k S (T i,j + (l/2) ) 



(B17) 


typify the notation for the conductivity terms in equations (B16) and (B19). 

The Tj j terms in equation (B16) must be linearized. Again, this is done by re- 
placing T? . by the first two terms in its Taylor expansion about - that is, by 

( T S “' 1, ) 4 + 4 ( t u 1) ) 3 ( t m- t m ") 


which reduces to 


4 ( t m' 1) ) 3t m - s ( T fi' 1> ) 

Substituting this last expression into equation (B16) and solving for T^ gives 

J 


T (m) = C l T i,j-l + C 2 T i+l, j 

U 


C 8 T ‘-AJ. + . C 4 ( T ?+1.J + + 6C 4 ( T 1 J~ 1) ) 4 * C 5 ,r i.Hl 

( T {“' 1> ) 3 + C 5 


C 1 + c 2 + c 3 + 8 C 


(B 18 ) 
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where 


C l- K i,j-(l/2) r j-(l/2)^T 

^ r B 

c 3 = k ( 

C 4= F i,i-l"( r i*J ar s)^ 

C 5 = \ i+ ,l/2)(^4 ir s)^ 


(B19) 


with A0 eliminated. 

Equation (B18) is in the point-Jacobi iterative form. Each point on the interface was 
relaxed by applying the successive overrelaxation scheme to the Jacobi method. 

Resolving the heat-transfer problem in the shield region involved the solution of a 
system of nonlinear two -point boundary value problems. Consider only the second equa- 
tion of equations (B7) since the following procedure is identical for the top shield (N ) ; the 
only difference is one of notation. Rearranging this equation (where i represented the 

1L 

1 n shield) gives 


W d_ 
rcr dr 



A i+1, r 


+ F. 


i,i-l 




+ F 




(B20) 


For simplicity, let 


G= W 
r a 


(B21) 


and 
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01 ( F i, i+l T i+l, r + F i,i-l T i-l, r 


(B22) 


Also for simplicity, the J verticals are renumbered so that J = 0 is the interface ver- 
tical and the last vertical is J = n (see fig. 32) . Substituting equations (B21) and (B22) 
into (B20) yields 



(B23) 


Equation (B23) was solved on each shield subject to the following boundary conditions. On 
the right (from eq. (Bll)), 


| = R digt = 0 (B24) 

r 

The left boundary of the shield region is the interface at r = R. . where the U values 

th lilt j_i 

are known at the m tn iteration as a function of the temperatures T at the m tn iteration; 
so on the left 



U ( R int,z) =f ( T > z ) (B25) 

The solution was effected by applying Kalaba’s linearization scheme (ref. 5) to 
equation (B23) and then integrating the result. The process just described was carried 
out as follows considering a general shield (i) . Using Kalaba’s method for linearization, 

L(U) = N(U) 


where L is a linear function and N is a nonlinear function, yields 

l(u^ = ^ 

i.1. it. 

where k represents the k tn iterate on the i cn shield. 

Regrouping superscripts (k) in equation (B26) gives 



(B26) 


(B27) 


where the LHS is a function of coefficients evaluated at the previous iterate and the RHS 
is also a function of the previous iterate. Applying equation (B27) to equation (B23) gives 
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L(U) 


= G— (r — \ + a? 
dr \ dr / 

N[T(U)] = + 

ff=<\l +1 + P M -l)Kr>] 


3 fLI 

dU 


where 


dT _ 1 

dU 


t s[' r < U i, r>] 


which follows from equation (B6) . Let 


■«-gr 


and 


b(x) 


= N (u (k_1) ) - 


dN^-^ydC-l) 
dU/ 


then, dropping superscripts k, equation (B23) becomes 


G— [ r — ) + a - a(x)U = b(x) 


dr \ dr 


(B28) 


I »Nl 


I = l 
1 = 0 


-Interface 


/ 


X- 

A 

Al,\) 

F * 

r 



J-0 J = 1 

Figure 32. - Renumbering of shield verticals. 


J = n 
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Consider the point (i, j) = (1, 1) where j s J (fig. 32). Integrating equation (B28) 
from j - ^ to j + ^ gives 


J ri+(l/2) r 

I -j'S**/ 

3 ~(l/2) Jj 


j+(l/2) /*j+(l/2) p+(l/2) 

a dr - I a(x)U dr = I b(x)dr 

-(1/2) Jj-(l/2) */j-(l/2) 


(B29) 


and using central differences for the derivatives reduces equation (B29) to 


r j ~(l/2) U j-1 


2r, 


a(Xj)(Ar g )' 


U 3 + P j+(l/2) U j+1 


U = [b( Xj ) - aj( Ar s ) ! 


G J 


(B30) 


Repeating the integration on equation (B28) for each point j, 1 < j ^ n, and using the 
boundary values at the end points result in a system of n equations in n unknowns. 

From equation (B30) and the statement following equation (B27), it can be seen that 

fkl 

this is a system linear in U v ' which can be written in matrix notation as 

AU< k ) = d 


where 


^ =^(r. ; 

is a square tridiagonal matrix of order n and 

d = d ^r. ; U^" 1 ^ 

is a column vector of n components . 

This system was also solved by the overrelaxation scheme applied to the line-by-line 
method. Each line (shield) was iterated on until the Euclidean norm of the U vector con- 
verged to within a tolerance of 0.01 percent. 

After each line in the shield region - beginning with the shield closest to the tank wall 
and moving outward to the uppermost shield - was converged, the next overall (m + 1) 
iteration was begun in the buffer region. This overall iteration was carried out until the 
relative residual error 
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I 


was less than 0. 001. Here 


|heat in | - |hea t out| 
|heat in 


heat in 


,N a ( T rs 

source 

boundary 


f * 

# rs 
J Radiating 



dA 



where E, F, and G are points on the tank surface shown in figure 32. 

In the iterative methods described, the overrelaxation parameter co is free to be 
chosen subject to the restriction 0 < o> < 2. It has been shown (ref. 4, ch. 4) for sys- 
tems of linear equations that there is an optimum value of o> for which convergence is 
fastest. Because of the nonlinearity of the problem being considered here, no rigorous 
analysis exists for the determination of a corresponding optimum value of co. Numerical 
experiments carried out at the Lewis Research Center, however, indicated the validity of 
the concept of an optimum value, and values of co obtained in these experiments were 
used in the production runs. The result for a typical 40-shield case was a> = 1.85 for 
1840 points leading to convergence of the relative error in the heat balance to less than 
0. 1 percent in 248 iterations. 
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